Long time deviation from exponential decay: non-integral power laws. 
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Quantal systems are predicted to show a change-over from exponential decay to power law de- 
cay at very long times. Although most theoretical studies predict integer power-law exponents, 
recent measurements by Rothe et al. of decay luminescence of organic molecules in solution {Phys. 
Rev. Lett. 96 (2006) 163601} found non-integer exponents in most cases. We propose a physical 
mechanism, within the realm of scattering from potentials with long tails, which produces a contin- 
uous range of power law exponents. In the tractable case of the repulsive inverse square potential, 
we demonstrate a simple relation between the strength of the long range tail and the power law 
exponent. This system is amenable to experimental scrutiny. 

PACS numbers: 



I. INTRODUCTION 

Rothe et al. [l[ have recently presented experimental 
evidence for "Violation of the Exponential Decay-Law at 
Long Times" . Deviations from exponential decay of an 
unstable quantum system had been predicted long ago as 
a result of the semiboundedness of the energy spectrum 
0. However, the transition into the non-exponential 
regime at long times occurs in most decaying systems 
after many lifetimes, where the minute survival proba- 
bility had previously frustrated verification. In addition, 
it had been argued that interaction with the environment 
and/or the measurement process would suppress the long 
time deviation [3| 1 . Rothe et al. measured the decay of 
luminescence of various dissolved organic molecules after 
pulsed laser excitation. The results clearly show a range 
of times over which the decay of the survival probabil- 
ity P(t) is exponential, followed by another where it is 
algebraic, P(t) ~ i - * 4 . 

One puzzling aspect of the experimental results of 
Rothe et al. [l[ is that their inferred algebraic expo- 
nents /j, are patently non-integral, except in one case 
(PtOEP, frozen solution), for which /i = 4. Other val- 
ues range from /j, = 2.08 to 4.07. These non-integer ex- 
ponents are in sharp contrast to theoretical studies of 
long-time deviation from exponential decay which predict 
or postulate integral exponents for the algebraic decay 
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1 In fact the perturbing effect of measurement depends on the par- 
ticular process. It has been shown to be negligible for distant 
detectors in solvable models 
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Reacting to this apparent contradiction, we explore 
physical mechanisms which could lead to long-time devi- 
ations from exponential decay with arbitrary, non-integer 
power law exponents. On the mathematical side, an "ex- 
planation" of an arbitrary power-law exponent follows 
from the expression for the survival amplitude in terms 
of an energy density oj{E). As defined in Section II, 
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where we have assumed no bound states and an energy 
threshold at E m = 0. The large-time asymptotics depend 
on the behaviour of uj(E) near the origin [l2j]. In fact, 
under mild conditions on u>(E) as a function of the com- 
plex variable E, A(t) ~ H"* 1 ) as t oo if uj(E) ~ E v 
as E -► 0. This implies P(t) = \A{t)\ 2 - t~ 2 ^ +1 \ 
From this perspective, non-integer values of the expo- 
nent fi = 2(y + 1) are the consequence of an appropriate 
non-integer v in the threshold density of states. 

By itself, this formal argument is not satisfying since 
it tells us nothing about the physical mechanism behind 
the threshold behaviour of w{E). We would like to learn, 
from the properties of the Hamiltonian or possibly of 
the initial state, the origin of such exponents, or, equiv- 
alently, the threshold behaviour of the energy distribu- 
tion. At this point we refrain from attempting a very 
general approach valid for the plethora of unstable sys- 
tems with exponential decay at intermediate times. In 
particular, the process of excitonic decay in dissolved or- 
ganic molecules is poorly understood, and present models 
explain only qualitati vely the processes that lead to their 
emission spectra (2ll . |22|, |23| . Hence we take refuge in a 
physical domain where full or sufficiently accurate ana- 
lytical and numerical methods are available, namely, the 
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decay of a wavepacket in a spherically symmetric poten- 
tial well, in a given partial wave. This is done with the 
hope that the lessons learned here may be useful else- 
where, although there is no guarantee that our results 
can be extended to such diverse fields as auto-ionizing 
states, photo-ionization, particle physics or spontaneous 
photon emission, because of the particular constructions 
of their distributions oj(E). 

One of the advantages of the decay of unstable states in 
ID and spherically symmetric 3D potential models is that 
it is a well studied subject, reviewed by Fonda, Ghirardi 
and Rimini some time ago [3j. The conditions leading 
to integer exponents are well understood. If the initial 
unstable state has angular momentum I and the poten- 
tial is "well behaved" at short and large distances, then 
the algebraic exponent is generally fi = 21 + 3. In par- 
ticular, for the most studied case of s-waves, fj, = 3; the 
i = contribution would be in any case the dominant 
one at very large times since it is the slowest decaying 
term. This algebraic decay law also applies in ID decay- 
ing systems |l4j. The derivation of these integral power 
laws will facilitate identification of the complementary 
cases in which the exponents may be non-integers, so we 
shall devote Section II to a review of the basic elements 
of potential scattering theory relevant for our purposes. 
(This presentation differs from Q in some respects, but 
the results are equivalent.) We shall conclude that poten- 
tials with long tails are natural candidates for providing 
non-integral decay exponents. The case of inverse square 
potentials, V(r) — 2m{3(0 + \)/(Ti 2 r 2 ) for r > r d , and 
[3 > (repulsive tail), is particularly suitable for theo- 
retical treatment and provides as a first approximation a 
very simple relation between the strength of the interac- 
tion and the exponent, /i = 2/3 + 3 (The attractive tail 
case (—1/2 < (3 < 0) will also be studied but the decay 
law is more involved.) The predicted algebraic exponent 
varies continuously with [3, with /i = 3 when (3 = cor- 
responding to a short range potential. This is shown in 
detail in Section III, for a potential model with inverse 
square tail in which to(E) is obtained analytically, and 
the conditions for the validity of the simple prediction 
for the algebraic exponent can be easily studied. In the 
final discussion we comment on the experimental feasi- 
bility of inverse square interactions. 



II. POWER LAWS FOR LONG TIME DECAY IN 
POTENTIAL SCATTERING 

We consider a system that is initially in a normalized 
non-stationary state | ^o) • The survival amplitude of that 
state is defined as the overlap of the initial state with the 
state at time t, and is the expectation value of the time 
evolution operator 

A(t) = (*o|exp(-*£rt/ft)|*o>, (2) 



H being the Hamiltonian. The survival probability, 
sometimes called "non-decay probability" [|[ , is 

P(t) = \A(t)\ 2 . (3) 

This is the function whose asymptotic time behaviour 
we will study. The stationary states of the Hamiltonian, 
H\®e,\) = E\§e,\), where by A we indicate any other 
quantum numbers characterizing states of energy E, de- 
termine a basis. In our models we have only a continuous 
spectrum, so the range of energies will be E m < E < oo, 
and the origin of energies is taken at threshold E m = 0. 
The completeness relation is then 
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which gives 
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cu(E) = ^K$ £ , A |* )| 2 . (5) 

A 

We will call u)(E) the "energy density" or "energy dis- 
tribution" of the initial state. To construct oj(E), we 
therefore need the continuum solutions denoted \<&e,\)- 
For a spherically symmetric potential, the Schrodinger 
equation has separable solutions $(f) = ug(r)Yf l (Q)/r. 
For each partial wave 

d 2 u e £(£+ 1) o 

TT- 2 ~ ue + [k 2 -v(r))u t = 0, (6) 
ar z r z 

where v(r) = (2m/h 2 )V(r) and k 2 = (2m/h 2 )E. For 
convenience, in the following we work with solutions 
wi(k,r) normalized as 

/•oo 

/ drwt(k',r)*w e (k,r) = 6(k' - k) (7) 
Jo 

and obeying the boundary condition 

lim w t (k, r) = (2/tt) 1/2 sin(fcr - nl/2 + SA (8) 

r — >oo 

where 8g is the phase shift for the partial wave £. These 
solutions are related to the regular solutions tfii (defined 
by their behaviour as Riccati-Bessel functions ji, 4>i(r) ~ 
ji(kr) when r — > 0) by 

\fiW\ 

where ft(k) is the Jost function. Among the different 
definitions we use the one by Taylor [24[ so that 

1 r°° 

Mk) = l+kJ drh t( kr H r )Mk,r), (10) 

see the Appendix A for a minimal account of Riccati- 
Bessel functions. 
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We now sketch the arguments that justify integral 
power laws for the long time deviation from exponential 
decay in scattering from "well-behaved" potentials, de- 
fined as those falling off faster than r~ 3 at infinity and by 
being less singular than r~ 3 / 2 at the origin [2J]. We shall 
assume for simplicity that the initial non-stationary state 
is in a particular partial wave £ without bound states, and 
localized, 

Ui(r) = for r > r a . (11) 
Then the survival amplitude takes the form 
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The function (fT0|) for £ > vanishes at k = if and 
only if there is a zero energy bound state whereas for 
£ = such an occurrence represents a zero-energy reso- 
nance [lij]. We disregard these exceptional possibilities 
and concentrate on the generic case fe{0) ^ 0. 

The main properties of concern for the long-time be- 
haviour of the survival probability are analyticity of <fie in 
the complex p-plane, its asymptotic behaviour near the 
origin (the same as the Ricatti-Bessel function) and at 
infinity [H], 

|^(fc^)l<7,( T ^)' +1 e |Im ^ 1 , (13) 

where "/i is some constant, and the behaviour of the Jost 
function at k = and at infinity. Eq. (fT0|) can be con- 
tinued analytically in the upper fc-plane (zeroes on the 
positive imaginary axis would represent bound states), 
and tends to one uniformly as \k\ — > oo in Im(fc) > 0. 

Now we deform the energy contour from the positive 
axis in Eq. (Q~2]) to one directed upwards along the lower 
imaginary axis on the first energy sheet. A long clos- 
ing arc at infinity does not contribute because the expo- 
nential growth of (/)£ implied in Eq. ([IB")) is limited by 
the localization (TlTj) and thus compensated by the de- 
caying time exponent. The behaviour of the integrand 
near threshold is thus, from Eq. (TTB"]) and the zero en- 
ergy asymptotics of the Jost function, of order E e+1 / 2 , 
so the long-time asymptotic of the survival amplitude is, 
according to Watson's lemma, A ~ t~(^+ 3 / 2 ) : and there- 
fore 

P(t) ~ r( 2i+3 \ (14) 

This is easily generalized to cases with a zero energy reso- 
nance or bound states, since the zero ener gy b ehaviour of 
the Jost function is known in these cases [24|, and again 
provides integer exponents. Another type of exception 
arises from the possible cancellation of the threshold de- 
pendence in Eq. (|13p by integration over the coordinate 
in (<j)\ui), by a suitable choice of uAf). In the context of 
ID potential models, Miyamoto has studied the long 



time behaviour of wavepackets scattered by a finite range 
potential. He showed that by careful adjustment of the 
low momentum components of the initial wave packet, 
initially located outside the interaction region, an expo- 
nent 5 for the decay power law may occur, instead of the 
generic value 3. As presented, such state manipulation 
produces only odd integer values. 

Comparing the result (fT4")) with the dependence on £ of 
the centrifugal barrier, it is tempting to speculate that a 
v(r) = /?(/? + l)/r 2 potential tail, where (5 > 0, may pro- 
duce long-time decay power laws of the form i~ ( 2 ^+ 3 ) for 
I = 0, based on the fact that the long tail plays for large 
r the role of a centrifugal term with a non-integer, effec- 
tive I. That this is indeed the case, will be demonstrated 
in a solvable model below. A proof of generality of this 
result is however far from trivial, since the properties of 
the functions used to arrive at Eq. (fl"4")) depend on fast 
fall-off conditions which are not satisfied for an inverse 
square tail. Instead of using £ = 0, we could alternatively 
assume from the start a non-integer, effective £ = (3 and 
its corresponding "partial wave" stationary equation to 
extend the usual results, but this requires a redefinition 
of the potential V(r) in the internal region to compensate 
exactly for the unphysical inverse square term. The con- 
sequence is an inverse square divergence in V(r) at r = 
which, again, makes the formal treatment non-standard. 
We shall therefore be content with a demonstration, via 
exactly soluble models, that the proposed long-time de- 
cay behaviour with non-integral exponents does indeed 
occur, and leave a more general and abstract theory for 
a further publication. 

III. DECAY OF AN UNSTABLE STATE IN 
POTENTIAL MODELS 

A. General results 

Consider a point particle of mass m in 3D subject to 
a central force. First we will derive some general results 
assuming only that beyond a certain radius, r = r^, the 
potential is inversely proportional to the square of the 
distance. We call this a "boundary condition model". 
For convenience, we write 

The inner part is also assumed to be free of singular- 
ities, even at the origin, so that the existence of regular 
solutions of the Schrodinger equation is guaranteed. To 
have a vanishing survival probability when t — > oo, we 
choose a potential with no bound states. Later, in order 
to have explicit analytic expressions for the properties 
of u>(E) and the asymptotic survival probability, we will 
use the more detailed "WB" (well-barrier) model. In that 
case, the inner part of V(r) will be a square well, depth 
— Vq < and extending from to r a , followed by a square 
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FIG. 1: Potential V(r) for one illustrative example: h = 2m 
1, vo = 0.5, v b — 1.8, r a — 3.0, r d = 3.4 and (3 — 1. 



barrier of height > located between r a and r&. An 
example is shown in figure [1] With suitable parameters 
for the well and barrier this model has well defined qua- 
sibound states, which lead to the exponentially decaying 
part of P(t). 

As announced in the Introduction, we will show that 
the value of j3 determines the exponent of the algebraic 
decay at long times. In our examples (3 > —1/2 will be 
either positive (repulsive), or negative but such that the 
attraction does not support bound state solutions. Since 
V(r) — > when r — > oo, the energy threshold in the 
model is indeed at E m = 0. 

For simplicity we take the initial state to have £ = 
0, and drop the subindex £ in the following. It is also 
convenient to define a different regular solution <p(k; r) — 
<p{k,r)/k that satisfies the boundary conditions 



0(fc;O) = O , 0'(fc;O) = l 



(16) 



and, as shown in provided that the potential goes 

to zero faster than 1/r when r — > oo, its asymptotic form 
is, see Eqs. (|515|) . 
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\f{k)\ 
k 



sin(/cr + 5). 



(17) 



where f(k) = \f(k)\e = /*(— k) on the real fc-axis. 
In addition, to be consistent with Eq. (U|), we define 
solutions normalized with respect to energy, we{t), 



w E (r) ■■ 
such that 
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n* k w(k '' r) ~-y n% 2 \f(k)\ 
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0(fc;r), (18) 



(w E \w E >) = S(E-E') = -^-6{k-k'). (19) 

h k 

Finally, we assume that the chosen potential allows a 
regular solution at zero energy, denoted 4>o{r) = 0(0; r). 



For the potential in Eq. (|15|) . the regular solution for 
r > can be written as a linear combination of Riccati- 
Bessel functions, 



(j){k,r) = am jp{kr) + b nl hpikr). 



(20) 



(The reason for the subindex /// is the consideration of 
two spatial regions for the inner part of the potential in 
the following subsection.) 

For repulsive potentials (3 > 0. We will also consider 
weakly attractive potentials with f3 > —1/2 so that the 
order of the related cylinder functions v — {3 + 1/2 re- 
mains positive; see Appendix A. For a detailed and up- 
dated discussion of the problems arising when the 1 /r 2 
potential becomes too strongly attractive, (3 < —1/2, see 
reference 28(. The constants am, bin, can be written 



in terms of the regular solution and its derivative at by 
matching the inner and outer solutions at this boundary, 

4>{k,r d ) = am jp(krd) + hii nfj(kr d ) 
dcj)(k, r) 



Or 



= k 



a/// j'p{krd) +bm n' {kr d ) ,(21) 



where the prime means derivative with respect to the 
argument of the Bessel function, x = kr. 

Using their asymptotic forms and comparing to Eq. 
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with 



C 2 = 
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J in 



= i(n' P ) 2 + (j' p ) 2 m,r d ) 2 + -^(k,r d ) 2 (n 2 +j 2 ) 



~<t>{k,rd)4>' {k,r d ){hph'p +j>%), 



(23) 



where, to simplify notation, <f)'(k,r d ) — (d(f>(r, k)/dr) rd , 
and the argument of the Bessel functions is x d — kr d . 
In subsection C we will need the asymptotic forms when 
k — > 0. In this limit, 4>(k, r d ) — > 4>o(r d ), which will always 
be finite and non vanishing unless we have made a rather 
unfortunate choice for the point r d . The value of the 
derivative will also be finite. Eqs. (f2"2")) and (p?3"|) then 
show that in this limit the explicit k dependence of the 
Jost function is determined by that of the Ricatti-Bessel 
functions: Eq. (fM)) gives an example of this. 

Since our initial state has I = 0, the sum over A in Eq. 
([5]) may be dropped, and writing the radial part of the 
initial state as \ui), the energy density is 



u(E) = \(w E \ 



(24) 



Once \ui) is chosen we can expand it in terms of the vari- 
ous solutions of the Schrodinger equation just described. 
We choose an initial state whose wavefunction is non- 
vanishing only when r < r a , 



Ui(r) 



sin k a r Q(r a - r), 



(25) 
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with k a = n a ir/r a and n a integer. In this paper we take 
n a = 1. Then ju,) coincides with the ground state of 
the inner well in the limit of infinite barrier height. This 
initial state is widely used in simulations of unstable sys- 
tems. Its advantage is its simple analytic form. At in- 
termediate times only the lowest quasibound state con- 
tributes significantly, which leads to exponential decay of 
P(t). For this choice of \ui) 
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which, combined with Eq. 
of the regular solutions. 



4>{k] r) sin k a r dr, (26) 



expresses u)(E) in terms 



B. The well-barrier (WB) model 

We now choose a specific model for the inner potential, 
consisting of an inner square well enclosed by a square 
barrier. In this case, 
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<p(k;r) = — smkjr , r < r a 
ki 

</}(k;r) = a n e Kr + b u eT Kr , r £ (r a ,r d ), (27) 

where kj — k 2 + vq, vo = 2mVo/h 2 , n 2 = Vb — k 2 and Vb — 
2mVb/h 2 . The constants an and bu are determined by 
matching at r = r a and from them we get the expressions 
for 4>{k; r d ) and its derivative, 



<!>{k;r d ) 



— | cosh 



ki 



ki 

KTb sin kjr a H sinh nrb cos kjr a I , 

K 

d>'(k;rd) — — ( sinh nr^ sin kjr a H — - cosh urb cos kjr a I , 
ki \ k 



with n =r d -r a . 



(28) 



Inserting these into Eqs. (f2"2"| and ([2^)) determines the 
Jost function. The overlap is now given by 



(w E \ui) 
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dr sin kjr sin k a r 
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and therefore when k a = n/r a (n a = 1), 



u{E) 



Airmk sin kir a 



h 2 rlk 2 \f{k)\ 2 [{*/r a ) 2 -k 2 ] 2 



(29) 



(30) 



Since k 2 



k' 2 



wo, even in this simple model the k- 
dependence is non-trivial. To illustrate our results, we 
fix on specific parameter values for our inner potential. 
We use units such that Ti 2 = 2m = 1 and set r a = 3.0, 
r d = 3.4, v — 0.5 and Vb — 1.8. Fig. [2] shows the 
computed energy spectra when the outer part of the po- 
tential has j3 ranging from —0.4 to 0.7. Fig. [3] shows the 




FIG. 2: Energy densities ui(E) v.s. energy for /3 = —0.4, 
—0.1,0.3 and 0.7, showing that peak heights increase with /3 
but lo(E) becomes singular for (3 < 0. 
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FIG. 3: Survival probabilities, P(t), exact (solid lines) and 
asymptotic approximation (dashed lines) of Eq. (|40p . From 
top to bottom: = -0.1, 0.3 and 0.7. 



corresponding survival probabilities. It can be seen that 
for this choice of inner parameters a strong resonance 
peak dominates the energy density in all cases, and that 
there is exponential decay of P(t) at intermediate times. 
Other smaller peaks present in lu(E) at higher energies, 
not shown in the figure, affect the corresponding P(t) 
only at very small times. We are not interested in that 
regime here. At long times the survival probability shows 
algebraic decay in all cases. To determine the form of 
P{t) when t —y oo, we need the asymptotic properties of 
w{E) near threshold. We find them now in the context 
of the boundary condition model. 



C. The threshold energy density 
From Eq. (JTHJ> , 

1 



w E {r) 
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(31) 
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and therefore 
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4>o(r) sin k a r dr. (32) 



Clearly in this limit all the ^-dependence is in the factor 
multiplying the integral, while the latter is to this order 
a constant. For the specific case of the WB model, the 
analytic expression for the integral can be obtained from 
Eq. (gSJ), replacing ki by fc/.o = y/vo- 

We look first at the behaviour of the Jost function, Eqs. 
(f2"2")l and (f2"3")l . also to lowest order in k. Using standard 
expansions, (see Appendix A,) the behaviour when f3 > 
as x = krd — ► is 



r l/2_ T /3+l 



2< 3 + 1 (/3+l/2)!' 
09-1/2)! 2^ 



7r l/2 a ,/3 



(33) 



and the contribution from the Neumann function will 
dominate in Eq. (f23|) . One easily finds that 



\f(k)\ ~ D fe-^ with 
p ^ 2^^- 1/2)! 



(34) 



where 0o = <f>o(i"d) and (/> = ^(r^). For the WB model, 
these are given explicitly in Eqs. (f2"5| . 
Inserting Eq. l|3"4")) into Eq. (j2"9"|) we arrive finally at the 
desired result, 



4mn sin 2 (fc/ j of a ) 
tfrl (kl-kjjkh P 2 



fc 2/3+l 



C fc 2/3+1 . 

(35) 



D. Decay at long times 



Usually, when an explicit expression for the energy 
density is available, to find analytic approximations for 
the asymptotic contribution to P(t) one replaces the in- 
tegration path in Eq. [5] by an equivalent one in the com- 
plex plane. Examples of these methods can be found in 
[12l |. [H and [T^. In our case, we choose a closed path 
consisting of the positive real energy axis, a quarter cir- 
cle of infinite radius joining the positive real axis to the 
negative imaginary axis, and the latter. As shown in the 
Appendix B, integration along the circular arc gives a 
vanishing contribution. Therefore 



poo 

A(t) = dE u(E) e~ lEt/h 
Jo 

/p—ioo 
dE oj(E) e ~ lEt/h + J dE uj(E) e 



= A p {t)+A v {t). 



iEt/h 



(36) 



We have analytically continued u>(E) into the lower half 
plane, and placed the cut running from the origin in the 



upper half plane. The term A p (t) has contributions from 
the poles enclosed in the contour. These are related to 
the resonances and lead to exponentially decaying terms, 
so we will not discuss A p (t) further. 

The second term, A v (t), gives the dominant contribu- 
tion to the decay at long time. To evaluate it, we make 
the change of variable E — —ix, 



A v (t) 



dx e xt ui(—ix). 



(37) 



In our examples the exponential decay becomes negligible 
beyond t ~ 200, so (assuming uj(—ix) is smooth), the 
range of values of x giving significant contributions to 
the integral is < x < x c ~ 1/200. Therefore, in a first 
approximation we use the small k expansion of u)(E), Eq. 
(f35|). to extend this function to the relevant part of the 
negative imaginary axis, and write 



(-0 



/3+1/2 C x /3+1/2 



(38) 



Inserting this into Eq. (|3T[) one immediately finds 

A v (t) ~ -H)^ +3/2 C T[fi + 3/2) i-W+3/ 2 ) (39) 
and therefore asymptotically 

P(t) ~ \A v (t)\ 2 = C 2 T 2 (/3 + 3/2) t- {20+3) . (40) 




FIG. 4: Effective /i/ v.s. (3. Continuous line with filled 
circles: standard choice of parameters for the inner part of 
the potential. Dotted line with filled diamonds: standard set 
except Vb = 1.6, dashed line with crosses: standard set except 
Td = 3.5. The straight dashed line is fj, = 2/3 + 3. 

This is the result promised in the Introduction: it 
shows that in this approximation the algebraic exponent 
is indeed fi — 2(3 + 3. Fig. [3] confirms graphically the 
accuracy of this result for a range of values of (3 either 
repulsive or weakly attractive. Note that this result is re- 
markable in two respects: because a) no other parameter 
of the potential affects this prediction; and b) the con- 
tribution of the higher order terms neglected in writing 
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Eq. (|33[) is negligible. Within the WB model we have 
checked that other choices of parameters for the inner 
part of the potential lead to the same conclusions. This 
will be further confirmed in figure 3J 

However, the above approximation fails for values of (3 
that make the outer potential attractive, especially when 
approaching the lowest value we considered, (3 = —1/2. 
This is shown in Fig. [5] for (3 = —0.4. Empirically we 
have found that for the range of times shown in the figure 
the asymptotic decay can still be well fitted by an alge- 
braic form, but the exponent determined from such a fit 
deviates from the simple law that holds for positive f3. 
Intuitively, the simple algebraic behaviour of the repul- 
sive case (3 > may be understood as a consequence of 
the dominance of the long centrifugal- like tail of the po- 
tential in near threshold scattering. However if the tail is 
atractive the particle is drawn to explore the inner part 
of the potential. This makes the decay more complex 
and also more sensitive to this inner part. We will now 
discuss that in further detail. 



with 
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Therefore 

uo{E) 



C k 2 



1 + (A /A_)fc 2 - + (A+/A_)fc 4 - 



(42) 



(43) 



replaces Eq. ([33)) . To have a result for A v {t) similar in 
form to the earlier one, we expand the denominator and 
write 




(44) 



600 



FIG. 5: Survival probability, P(t), for (3 = -0.4, (u = 0.1). 
Continuous line: exact; dashed line: single term asymptotic 
approximation Eq. (|40|l : crosses: the expansion of Eq. (I45p 
truncated at 4 terms. 



Looking again at the expansions of the Ricatti-Bessel 
functions in Eq. (|A2[) one sees that when j3 — > —1/2 the 
series for jp and hp both begin with x 1 / 2 and therefore 
the contribution from the jp cannot be discarded. At 
this stage it is convenient to use the auxiliary index, v — 
(3 + 1/2, that corresponds to the index of the cylindrical 
Neumann and Bessel functions in Eq. (|A1[) . We will 
be interested in small values of v. Our representative 
example will be that of Fig. [5l v — 0.1. Including only 
the lowest order terms, p = in Eq. (|A2|) . but for both 
functions, one finds 



m— 1 



with £2 = C already defined, and the other C^m obtained 
from Eqs. (|42|) and f|43[) . Repeating the derivation in 
Eqs. (J3HD to (gDJl, we obtain 



A v (t) ~ - J2 C2m r(l + mv) (it)-** 1 



-\-mv) 



(45) 



\f{k)\ 2 ~k(\_k- 2v + \ + \ + k 2v ) 



(41) 



The first term in this expansion reproduces Eq. 
but when v ~ 0.1 the algebraic exponents of the first few 
terms in the series will be comparable, and keeping only 
the first would be insufficient. Confirming this, figure [5] 
shows that with a single term in the series the prediction 
is wrong by an order of magnitude, and that one has to 
include at least four terms in the series to reproduce the 
exact P(t) for t > 200. When we use the approximation 
for uo{E) of Eq. (T4"3"|) in Eq. (|3"T|) , the agreement is even 
better. 

These results show that for positive (3 one can neglect 
the ^-dependence in <p(k,rd) and its derivative, and also 
the higher order terms in the series for the Ricatti-Bessel 
functions, and still have a satisfactory approximation for 
the asymptotic decay of the survival probability over the 
range of time included in the figures. But when (3 be- 
comes negative (outer potential attractive) the trunca- 
tion leading to the asymptotic law in Eq. (|4"0"]) ceases to 
be valid. Even then, we have empirically found that the 
exact P(t) can be well fitted with a formally similar alge- 
braic expression in a range of times that we will specify 
as the interval 400 < t < 800. To be more precise, and to 
try to simulate the way in which the algebraic exponents 
could be extracted in experiments like that of Rothe et 
al, we have made linear least square fits to the exact val- 
ues of \nP(t) v.s. ln£ for equally spaced values of t in the 
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above range. The fits are always excellent, so that in that 
range of times one can hardly distinguish in the figures 
the exact P(t) from the fitted form P f (t) = M . The 
effective algebraic exponents, /i/, thus found are shown 
in figure [4j One sees that when 8 < the effective ex- 
ponents depend on the parameters of the potential, and 
deviate from the simple law [i = 28 + 3. Making the 
barrier wider or higher, increases the values of /if for a 
given P, but in all cases as f3 — > —1/2 one finds that the 
values of \if become constant. 

We tried to derive explicit analytic expressions for P(t) 
in the limit v = 0. Using only the first order term in the 
expansion of the Neumann function, Nq(z), the first term 
in the asymptotic series for A v (t) goes like l/[t ln(i/r^)] 
with the next order terms involving higher powers of the 
logarithm. However, using only the first order term of 
the series for Nq is not enough; one needs at least one 
more term in both Jo and Nq. With these added terms, 
we were unable to find a useful analytic expression for 
A v (t). Still, the fact that the lowest order results include 
hit and powers thereof indicates that the algebraic form 
of decay for the asymptotic part of P{t) is not as univer- 
sal as previously thought, ft is only for potentials with 
a vanishing or repulsive outer part that such a simple 
algebraic decay law can be established. 



deviation is altered but not in the sought for algebraic 
form. 

Despite our retreat from the arena of the organic 
molecule experiments, the long-time decay described here 
is not simply academic. Attractive inverse square po- 
tentials occur physically as effective radial potentials be- 
tween a charged wire and a polarizable neutral atom [30| , 
the strength factor being proportional to the square of 
the linear ch arg e density of the wire and thus control- 
lable [10, [H], H3 ■ Combined with a repulsive centrifugal 
term, an arbitrary a/r 2 potential may be implemented. 
In addition, it is possible to modify the inner region and 
implement a potential minimum by a time varying si- 
nusoidal voltage in the high frequency limit [32], or by 
replacing the wire by a charged optical fiber with blue de- 
tuned light propagating along the fiber and the cladding 
removed 31| . Decay experiments with cold atoms show- 
ing exponential laws have been performed [30j |. and the 
ability to modify the potential parameters makes the ob- 
servation and study of the long-time power-law in these 
systems a realistic prospect. 
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Our study of the long-time deviations from exponen- 
tial decay, was motivated by the apparent contradiction 
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sults which unambiguously show such deviations, and 
theoretical models: the power laws in the experiments 
have non-integer exponents whereas the models predict 
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integral exponents? The complexity of the system stud- 
ied experimentally (large, excited organic molecules in 
solution), made us consider, instead of an ab initio or re- 
alistic approach, the more modest goal of understanding 
and answering the question in a tractable system. 

Limiting ourselves to a single-particle model of wave- 
packet decay from a scattering potential region, we have 
established a link between non-integer exponents and 
long potential tails, which is expressed by a simple for- 
mula for repulsive inverse square potentials. Before ar- 
riving at these results we explored several other possible 
sources of "anomalous" long-time decay: in particular, 
we tried to simulate the effects of the solvent molecules 
and of temperature by introducing randomly fluctuating 
perturbations in the potential and averaging over many 
realizations. This lead to localization effects and a non- 
vanishing constant value of P{t) at long times but non- 
integral exponents were not found. Another attempt was 
consideration of complex potentials that could simulate 
the effect of coupled channels due to measurement or 
other effects. However the result of an imaginary term is 
similar to that described in Eq. (11): the long-time 
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APPENDIX A 



The Riccati-Bessel functions jp(x),fi/3, and h^(x) — 

hp±ijp are defined as x times the corresponding spherical 
functions, or, in terms of cylindrical Bessel and Neumann 
functions, 



np{x) = J —N 0+1/2 {x). 



In this paper f3 may be non-integer. 
Their series representations are [29j . 
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U 2 J JPK ' sm(8+ 1/2)tt 
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APPENDIX B 

We will prove here that the energy density extended 
to the quarter circle at infinity, in the fourth quadrant of 
the energy plane, is vanishingly small and therefore does 
not contribute to A(t) in Eq. ([36]). To do so, we first 
rewrite the energy density as 



w{E) 



Arrnt 



h 2 ^rlk]({Tr/r a y-k}y 
sin 2 fc/r a 
C 2 ' 



G(E) 



(Bl) 



We will now extend these functions to the fourth quad- 
rant in the energy plane and look at their asymptotic 
behaviour when \E\ — > oo. It is convenient to write 



E = \E\e 2lv , 
k = |fc|e^, 



(B2) 



and ip £ (— 7r/4, 0). Similarly we will write fcj = \ki\e l>Pl 
and will now find an approximate expression for the latter 
valid in the limit when \k\ — » oo. From kj = k 2 + vq, we 
can write 



\ki\ 2 e 2i ^ = \k\ 2 e 2lv +v 



|fc|e l¥ 



*■'() 



|fc/|e**" 
ki 



2|fe|e*> 



vo_ * 
2\k\ 



vq 

2\k\' 



o(\k\- 2 ) 
o(\k\- 2 ). 



Similarly, since k 2 = k 2 — Vb, 



kh 



Vb 

2\k\ 



e-^ + 0(\k\- 2 ) 



(B3) 



(B4) 



In the same limit, the sine in the numerator of G(E) will 
give 



sin kjr a 



1 

27 f 



Ak 7 r c 



(B5) 



since Im(fc/) < 0. Let us now consider the denomina- 
tor: from Eq. (f23|) one sees that both the Ricatti-Bessel 
functions and the regular solutions at r = must be 
extended to complex energies in the fourth quadrant and 
in the large \E\ limit. Using the asymptotic expansions 
given e.g. in [331 ]. we find 



■n 



I \2 



1) 



{n'pY + Up) 



2z 2 
2z 3 



0{z~ 



with z — krd 
order z~ 4 



2z 2 

Inserting these into Eq 



+ 0(z- 5 ) 

[ l + 0(z-% (B6) 



we find to 



,i p{p_ 
k 



p(p + 1) 



2z 2 



1 



1) 



2z 3 



(B7) 



Now one has to deal with <f> and cf>' , whose explicit expres- 
sions are given in Eqs. 1281 Using the approximations in 
Eqs. (|B3l EU) and (|B5ll. one finds 



AM 2 



k 2 



1 

k 2 



1 

¥' 



_Vb_+VQ 2i(k b r b +k ira ) 

4fc 4 

_^_2i(k b r b +k ira ) 
fc 4 

_A_ e 2 i(k b r b +k ira )^ ^ 



so that all terms in the denominator of Q (E) diverge like 
exp(2i(fc/r a + fc(,r;,)). Since the sin 2 in the numerator 
diverges only like exp(2zfc/r a ), see Eq. (|B5I) . the result is 
that |£7(£0| — > and therefore the contribution to A(t) 
due to the integral on the quarter circle at infinity is nil. 

Note that this result can be easily extended to any 
value of n a introduced in Eq. (|25[) . Therefore any initial 
state, \ui), that can be expanded as a finite superposi- 
tion of eigenstates of the infinite square well will have an 
energy density with similar properties. 
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